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Dynamical systems often contain oscillatory forces or depend on periodic potentials. Time or space 
periodicity is reflected in the properties of these systems through a dependence on the parameters of 
their periodic terms. In this paper we provide a general theoretical framework for dealing with these 
kinds of systems, regardless of whether they are classical or quantum, stochastic or deterministic, 
dissipative or nondissipative, linear or nonlinear, etc. In particular, we are able to show that simple 
symmetry considerations determine, to a large extent, how their properties depend functionally on 
some of the parameters of the periodic terms. For the sake of illustration, we apply this formalism 
to find the functional dependence of the expectation value of the momentum of a Bose-Einstein 
condensate, described by the Gross-Pitaewskii equation, when it is exposed to a sawtooth potential 
whose amplitude is periodically modulated in time. We show that, by using this formalism, a small 
set of measurements is enough to obtain the functional form for a wide range of parameters. This 
can be very helpful when characterizing experimentally the response of systems for which performing 
measurements is costly or difficult. 

PACS numbers: 05.45.-a, 05.60.-k, 05.60.Cd 


It often happens that a system—physical or 
otherwise—can be described using a model that in¬ 
cludes one or several periodic functions, with the same 
or different periodicities. Whether these functions 
represent external oscillatory forces, modulating ampli¬ 
tudes, or space-periodic potentials is immaterial for our 
forthcoming discussion, as are the specific details of the 
underlying dynamics, which can be either deterministic 
or stochastic, classical or quantum, dissipative or nondis¬ 
sipative, linear or nonlinear, etc. For our purposes, the 
only relevant feature that all these systems have in 
common is that their properties depend in a certain way 
on the periodicities, amplitudes, and relative phases of 
these functions. 

Examples of these types of systems abound in the lit¬ 
erature. Because of their ability to describe a wide vari¬ 
ety of phenomena, the most significant ones are probably 
the periodically driven systems [T]. To name but a few 
instances, these driven systems have proven to be use¬ 
ful in the study of stochastic resonance [2], vibrational 
resonance [3], classical and quantum stochastic synchro¬ 
nization [3] , opinion formation processes , coherent de- 
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struction of tunneling , dynamical localization and de- 
localization [ 7 ], ratchet effect [8], and atomic quantum 
motors [S] . It is worth mentioning that some of the above 
examples (e.g., the ratchet effect) include both time- and 
space-periodic functions. 

The purpose of this paper is to provide a general de¬ 
scription of this widespread situation, namely, a system 
whose dynamics depends on a set of periodic functions, 
and to discuss the consequences of some symmetries that 
are often encountered in this class of systems. Specif¬ 
ically, we show that the functional dependence of the 
system’s properties on some of the parameters contained 
in these functions may be determined to a large extent 
by simple symmetry considerations. 

To be more precise, let us suppose that the dynamics 
of the system under consideration depends on N periodic 
functions of a single variable (j. This variable may repre¬ 
sent time, space, or even a generalized coordinate of some 
sort. As a matter of fact, our forthcoming discussion can 
be readily extended to the case in which there are more 
than one variable—say, space and time, or several spatial 
variables—, but for the sake of clarity we first consider 
only one. These periodic functions are assumed to be of 
the form 

fjiO = otj coslUljiC - Co)] + Pj sin[fIj(C - Co)] , (1) 

with j = I,..., A^, where Ulj are the (temporal, spatial, 
or generalized) angular frequencies, and aj and /3j the 



2 


partial amplitudes. This assumption is not as restric¬ 
tive as it might appear, since any well-behaved periodic 
function can be approximated to any desired degree of 
accuracy by a finite sum of trigonometric functions. The 
parameter Co has been introduced in Eq. Q to simulta¬ 
neously shift all the functions fj{C) along the C-axis. 

Let T represent a certain (physical) quantity of the sys¬ 
tem. We are interested in the functional dependence of 
T on some of the parameters appearing in the functions 
/j(C). To study this dependence, we will make use of 
the following simple rule: any transformation of the pa¬ 
rameters defining the functions fj{C)t which leaves these 
functions invariant, also leaves the value of the quantity 
T invariant. For this to be true, it is evidently assumed 
that all the other parameters in the problem remain un¬ 
changed during this transformation. 

In order to apply this rule, it is convenient to rewrite 
the periodic functions in Eq. Q in the form fj{C) = 
ej cos [fIj(C — Co) + ‘Pj], where we have introduced the N 
amplitudes Cj = (a^ + and the N phases (pj, 

satisfying the equations cospj = cxj/tj and sinc/j^ = 
—Yjltj. Let us define the vectors = (fli,..., IIat), 
if = {pi, ..., Pn)^ and e = (ci,..., e^r). Then, it is clear 
that the set of periodic functions is invariant under these 
two transformations 

Ti : {Co,rj,¥>,e} I—^ {(^o,n,p + TTU^^\e^^Y , (2) 

T 2 : {Co,^^,¥’,e} '—^ {0, rj, - CoJ^, 4 ) (3) 

where denotes the jth row of the N x N identity 
matrix and the vector is obtained from the vector e 
by replacing its jth component by —Cj. Consequently, 
from the above-mentioned rule, 

T(Co, P, e) = T(Co, P + = T(0, p - ^Co, e), 

(4) 

where we have only explicitly written the dependence of 
T on the parameters Co, P, and e, for fl is assumed to 
be fixed throughout this study. The analysis of the de¬ 
pendence on ft requires the use of alternative techniques 

By applying the first equality in the above equation 
twice, we see that T(Co,¥’,e) is periodic with respect 
to all the components of the vector p with period 27r. 
Therefore, taking into account the second equality in 
Eq. Q, it can be expanded in Fourier series as 

T{(:o,P,e)= (5) 

feez« 

where p ■ k = Pj^j and 

u44 =![■■■ J[ ^ Tf (0, . (6) 

Without loss of generality, we can assume that the 
quantity T is real. Otherwise, one would consider its 
real and imaginary parts separately. Then, taking into 
account that the imaginary part of Eq. ([^ is zero and 


introducing the functions 7 fe(e) = 'yfc(e)n^i S' 
obtains 

^iCo,P,e)= l7fc(e)|("]Jef"'^ 

\i=i / 

xcos[{p-flCo)-k + Xk{€)], ( 7 ) 

where | 7 fc(e)| and Xfc(e) are, respectively, the modulus 
and phase of the complex number 7 fe(e). Note that, ac¬ 
cording to Eq. § and the first equality in Eq. ffl, the 
functions 7 fc(e) and, consequently, | 7 fc(e)| and are 

even in each of the arguments ej. In addition, since T is 
real, | 7 fc(e)| = | 7 _fc(e)| and . 

An important result follows from assuming that the 
quantity T is invariant under arbitrary shifts of all the 
periodic functions fj{C) along the ^-axis. In this case, T 
is independent of Co and, accordingly, all the coefficients 
7 fc(e) such that k ■ ft Y ^ must be zero. Thus, dropping 
the dependence of T(Co, P, e) on Co, we obtain 

N 

'^{p,£)= Y l7fc(4|cos[v3-fe-kXfc(e)]n4^''> 

keSn i=l 

where Sn is the set of all ordered Wtuples of integers 
orthogonal to ft, i.e., iSfj := {fc G : ft ■ k = 0}. It 
is worth noting that, independently of whether or not T 
satisfies the above-mentioned shift-invariance property, 
Eq. (|^ is always valid for the average value T{p,e) = 

limACo-*-oo(ACo)"^ /(fdCo T(Co, P, e). 

Let us now consider the case in which the N angular 
frequencies fli,... ,fl]\r are incommensurable, i.e., it is 
not possible to express one of them as a linear combina¬ 
tion of the others with rational coefficients. Then, the set 
Sfi consists of the single element fe = 0 and, according 
to Eq. ([^, T{p,e) = 70 (e). Consequently, the quantity 
T is an even function in each of the components of the 
vector e and is independent of the phases p. If there 
were additional symmetry considerations leading to the 
conclusion that T is an odd function in any of the compo¬ 
nents of e, then necessarily T = 0 (see, e.g., Refs. fTTlIT^ 
for the case of two incommensurable frequencies). 

In contrast, if the N angular frequencies are pairwise 
commensurable, then there exists a frequency ft such that 
ft = fin, where n = (ni,..., un), with rij being positive 
integers. Thus, the condition ft ■ k = 0 is equivalent to 
the Diophantine equation n ■ k — 0. The general solu¬ 
tion of this last equation can be expressed as an integer 
linear combinations of a set of N — 1 generating vectors, 
k^^\ ..., k''^~^\ each of which satisfies n ■ k^^^ = 0 |I3) . 
Consequently, Eq. (|^ now reads 

T(¥>,e)= ^ l7fc(q)(4l 

qez«-i \i=i / 

X cos [p ■ fe(q) -k Xfc(q)(e)] , (9) 
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where fe(q) = Hence, T depends on 

only through the collective phases (p-k^^\ ... ,(p- . 

Equation ^ can be considered as a generalization of 
the results reported in Ref. m in the context of rocking 
ratchets. 

We now proceed to study the perturbative behavior of 
the quantity T for sufficiently small values of the ampli¬ 
tudes Cj, expressed in suitable dimensionless units. To 
this end, we will assume that, for Co = Oj the quan¬ 
tity T, expressed as a function of the partial ampli¬ 
tudes a and ,3, has a Taylor power expansion of the 
form ^ with No being the set 

of all nonnegative integers and ai^r the coefficients of 
the Taylor series. Then, using that ttj = ej cos (pj and 
Pj = —ejsiinpj, as well as Eq. and the definition of 
the functions 7 fc(e), we find 


N 


7fc(e) = 


lj+rj-\kj\ 


l,rGN" 


i=i 


( 10 ) 


where 


^Lr,k — 


— COS* ip sin’" ip e . 

27r ^ ^ 


( 11 ) 


The above integral vanishes unless I + r — |fc| is a non¬ 
negative even integer. Consequently, 


N 


7k(e)= 

pSN™ i=i 


( 12 ) 


with 


N 


bk.p — ^ 0,l^rY\_^h^rj,kjS2p^, 


lj+rj-\kj\ ■ 


(13) 




The series expansion in Eq. (12 1 , together with Eqs. 

([^, or (§, depending on the specific case, allows deter¬ 
mining the functional dependence of the quantity T on 
the parameters Coj ^md ip. In practice, for sufficiently 
small amplitudes, these expansions can be truncated to 
include only a few terms. In this case, the determina¬ 
tion of T(Co,¥’,e) is reduced to the calculation of a few 
model-dependent coefficients b^ p. 

In order to understand how to put these ideas into 
practice, let us consider a one-dimensional Bose-Einstein 
condensate described by the nonlinear Gross-Pitaevskii 
equation [15] 


ih 


d'S{x,t) 

Wt 


2m dx^ 


U{x,t)'i>{x,t) 


+ g\'i’{x,t)f-^{x,t ), 


(14) 


where 4'(a;,t) is the condensate wave function, normal¬ 
ized to 1 , TO the mass of the bosons, g the scaled strength 
of the nonlinear interaction, and U{x,t) a potential of 


the form U{x,t) = [/oH(a;) [1-I-cos(27rt/T)], with Uq 
being a constant with the dimensions of energy and 
T the temporal period. The spatial part V{x) is de¬ 
scribed by a biharmonic function of the form V(x) = 
cos(27ra;/A) -I- r? cos( 47 ra;/A -I- </>) [T^, where g is the am¬ 
plitude of the second harmonic relative to that of the 
fundamental, A the spatial period, and 0 the relative 
phase between the two harmonics. In order to numer¬ 
ically solve Eq. ( [l4| ), we consider an initial condition of 
the form 'I'(a:, 0) = \/2/Acos(27ra;/A), and impose peri¬ 
odic boundary conditions [T^, i.e., 4'(a; -b A,t) = 4'(a;,t) 
Vx € K. In this paper, we focus our attention on the 
expectation value of the momentum evaluated at time 
t = T, which is given by the expression 

Pt = J dx'i>*{x,T)(^-ih^^{x,T). (15) 


Specifically, we are interested in the functional depen¬ 
dence of Pt on the parameters </> and g, i.e., in the func¬ 
tion Pt(^, g)- 

To obtain an expression for PT{(j)^g), we first apply 
the previously developed formalism to the periodic func¬ 
tions fi{x) = cos(27ra;/A) and f 2 {x) = r 7 COs( 47 ra;/A + (j)). 
The correspondence with our previous notation is C = a:, 
(o = 0, = (27r/A,47r/A), if = (0,0), and e = ( 1 , 77 ). 

Consequently, using Eq. 0 together with Eq. ( [T^ leads 
to 


k,p—0 

(16) 

where fXk,p and k'k^p are coefficients that do not depend 
on 0 and g. Since Eq. ( jldj ) as well as the consid¬ 
ered initial condition are invariant under the transfor¬ 
mation Ta : {a:,0} '—> {—a;, —0}, it is easy to see that 
Pt(—0, 77 ) = —Pt( 0, ??)■ Hence, all the coefficients Vk.p 
in Eq. (16 1 must vanish. 


Let us now assume that g is sufficiently smal l so that 
we can neglect the terms of order 0 ( 77 ^) in Eq. (161, and 
approximate Pt(0, t?) by 

Pt(0, 77 ) « 77 (/ri,o + 77 ^ 171 , 1 ) sin 0 -b 77 ^/ 72.0 sin( 20 ). (17) 

The determination of the function Pt(0, 77 ) is thus re¬ 
duced to evaluating the three coefficients pig, Pi.i, and 
P 2 ,o- These coefficients can be easily calculated if we 
know, e.g., PT( 7 r/ 2 , 771 ), PT( 7 r/ 2 , 772 ), and PT( 7 r/ 4 , 771 ), 
with 771 and 772 being two different values of g within 
the validity range of Eq. (dzl)- In that case, 

77|PT(7r/2,77i) - 77fPT(7r/2,772) 


Mi.o — 


771772 ( 7 ?! - gj) 


(18) 


and 


b2,0 = 


7?2PT(7r/2,77i) - 77iPT(7r/2, 772) , . 

2 Pr( 7 r/ 4 , 771 ) - \/2gi {pifi + glpip) 


2gl 


( 20 ) 
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FIG. 1. (Color online) Dependence of the dimensionless mo¬ 
mentum Pt{4>,v) ~ XPt{4’,v) on fhe relative phase (f) for 

rj = 0.1, 0.3, 0.5, 0.7, and 0.9. The rest of the parameter 
values, expressed in dimensionless form, are h~^X^T~^m = 1, 
h~^TUo = 1, and h~^X~^Tg = 1. The results obtained from 
the numerical solution of the Gross-Pitaevskii equation in 
Eq. (141 are represented by five different types of symbols. 


as shown in the upper legend box. The three solid symbols 
indicate the values used to evaluate the coefficients pi.o, Pi,i, 
and g, 2 ,o from Eqs. ( |18[ ), ( |19[ ), and ( |20[ ). The results calculated 
by using Eq. (171 with the obtained coefficients are depicted 
with five different types of lines, as indicated in the lower 
legend box. 


ues of PTi4’,v) into Eqs. (18), (19), and (20), one can ex¬ 
plicitly calculate the coefficients that appear in Eq. 0. 
yielding H ^Xfii^o ~ 6.1571, h « —1.2083, and 

h~^X^ 2 ,o ~ —0.6417. The results for PT{(j),ri) obtained 
by using Eq. 0 and these three values of the coeffi¬ 
cients are depicted with five different types of lines in 
Eig.g To check the accuracy of these predictions, we 
also show in Eig. [T] with five different types of symbols, 
the results obtained from the numerical solution of the 
Gross-Pitaevskii equation in Eq. (14). As can be seen. 


the agreement between our predictions and the numerical 
results is excellent. 

In conclusion, we have developed a general theoretical 
framework for describing dynamical systems that contain 
periodic terms. The formalism can be equally applied 
whether the system is classical or quantum, stochastic or 
deterministic, dissipative or nondissipative, linear or non¬ 
linear, and more importantly, regardless of whether the 
periodic terms are time oscillations (in external forces or 
modulating amplitudes) or periodic spatial potentials (as 
in the example), or both. We have shown that the func¬ 
tional dependence of the system’s properties on some of 
the parameters of the periodic terms can be determined, 
to a large extent, by simple symmetry considerations. In 
particular, within the appropriate range of parameters, 
this functional dependence can be obtained, except for a 
few unknown constant coefficients. This last result can 
be very helpful when characterizing experimentally the 
response of systems for which performing measurements 
is costly or difficult. 


To sum up, it is enough to know Pt for three different 
values of (</>, rj) to determine its value for a wide range of 
parameters </> and rj. 


In order to illustrate this result, we have used a 
spectral method to numerically solve Eq. ( [l4| ) for 
(^, 77 ) = ( 77 / 2 , 0 . 1 ), ( 77 / 2 , 0 . 3 ), and ( 77 / 4 , 0 . 1 ). The rest 

of the parameter values, conveniently expressed in di¬ 
mensionless form, are h~^X'^T~^m = 1, h~^TUo = 1, 
and h~^X~^Tg = 1. The results obtained for the di¬ 
mensionless momentum PT{4>,rf) = h~^XPT{4>,r]) are 
PT(7r/2,0.1) « 0.6145, PT(7r/2,0.3) « 1.8143, and 
iV(7r/4, 0.1) « 0.4281. These three values are shown by 
solid symbols in Eig. By substituting these three val- 
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